An Intuitive Electric-field Contribution Decomposition Model for Chemical Processes and Its Applications on Diels-Alder Reactions

External electric field(EEF) has shown its advantages in tuning chemical reaction as an efficient and feasible-to-control tool. In this paper, we explored the mechanisms of three EEF-regulated Diels-Alder reactions including two traditional-DA reactions to form two C-C single bonds and a hetero-DA reaction to form both a C-C and a C-O bond, respectively, and introduced an EEF contribution decomposition(ECD) model to understand how the EEF coupled with the intrinsic nuclear and electronic redistributions so as to affect chemical reaction. The ECD model, by decomposing the overall EEF effects into geometry re-equilibrium and static induction parts, can give a clear and quantitative picture of a physical quantity change upon EEF, as demonstrated on relative energies, activation barriers, charge distribution and dipole moments. The ECD analyses will shed light on the effective tuning of chemical reactions by the electric field.     

Swendsen‐Wang dynamics for general graphs in the tree uniqueness region

The Swendsen‐Wang (SW) dynamics is a popular Markov chain for sampling from the Gibbs distribution for the ferromagnetic Ising model on a graph G = (V,E). The dynamics is conjectured to converge to equilibrium in O(|V|^1/4) steps at any (inverse) temperature β, yet there are few results providing o(|V|) upper bounds. We prove fast convergence of the SW dynamics on general graphs in the tree uniqueness region. In particular, when β < β_c(d) where β_c(d) denotes the uniqueness/nonuniqueness threshold on infinite d‐regular trees, we prove that the relaxation time (i.e., the inverse spectral gap) of the SW dynamics is Θ(1) on any graph of maximum degree d ≥ 3. Our proof utilizes a monotone version of the SW dynamics which only updates isolated vertices. We establish that this variant of the SW dynamics has mixing time and relaxation time Θ(1) on any graph of maximum degree d for all β < β_c(d). Our proof technology can be applied to general monotone Markov chains, including for example the heat‐bath block dynamics, for which we obtain new tight mixing time bounds.     

Jordan–Cattaneo waves: Analogues of compressible flow

We review work of Jordan on a hyperbolic variant of the Fisher–KPP equation, where a shock solution is found and the amplitude is calculated exactly. The Jordan procedure is extended to a hyperbolic variant of the Chafee–Infante equation. Extension of Jordan’s ideas to a model for traffic flow are also mentioned. We also examine a diffusive susceptible–infected (SI) model, and generalizations of diffusive Lotka–Volterra equations, including a Lotka–Volterra–Bass competition model with diffusion. For all cases we show how a Jordan–Cattaneo wave may be analysed and we indicate how to find the wavespeeds and the amplitudes. Finally we present details of a fully nonlinear analysis of acceleration waves in a Cattaneo–Christov poroacoustic model.     

More efficient time integration for Fourier pseudospectral DNS of incompressible turbulence

Time integration of Fourier pseudospectral DNS is usually performed using the classical fourth-order accurate Runge-Kutta method or other second- or third-order methods, with a fixed step size. We investigate the use of higher-order Runge-Kutta pairs and automatic step size control based on local error estimation. We find that the fifth-order accurate Runge-Kutta pair of Bogacki and Shampine gives much greater accuracy at a significantly reduced computational cost. Specifically, we demonstrate speedups of 2× to 10× for the same accuracy. Numerical tests (including the Taylor-Green vortex, Rayleigh-Taylor instability, and homogeneous isotropic turbulence) confirm the reliability and efficiency of the method. We also show that adaptive time stepping provides a significant computational advantage for some problems (like the development of a Rayleigh-Taylor instability) without compromising accuracy.     

Investigation on the effect of density ratio on the convergence behavior of partitioned method for fluid–structure interaction simulation

In order to investigate the effect of density ratio of fluid and solid on the convergence behavior of partitioned FSI algorithm, three strong-coupling partitioned algorithms (fixed-point method with a constant under-relaxation parameter, Aitken’s method and Quasi-Newton inverse least squares (QN-ILS) method) have been considered in the context of finite element method. We have employed the incompressible Navier–Stokes equations for a Newtonian fluid domain and the total Lagrangian formulation for a non-linear motion of solid domain. Linear-elastic (hyper-elastic) model has been employed for solid material with small (large) deformation. A pulsatile inlet-flow interacting with a 2D circular channel of linear-elastic material and a pressure wave propagation in a 3D flexible vessel have been simulated. Both linear-elastic and hyper-elastic (Mooney–Rivlin) models have been adopted for the 3D flexible vessel. From the present numerical experiments, we have found that QN-ILS outperforms the others leading to a robust convergence regardless of the density ratio for both linear-elastic and hyper-elastic models. On the other hand, the performances of the fixed-point method with a constant under-relaxation parameter and the Aitken’s method depend strongly on the density ratio, relaxation parameter selected for coupling iteration, and degree of deformation. Although the QN-ILS of this work is still slower than a monolithic method for serial computation, it has an advantage of easier parallelization due to the modularity of the partitioned FSI algorithm.     

Aeroelastic analysis of cantilever plates using Peters’ aerodynamic model,and the influence of choosing beam or plate theories as the structural model

In this paper, the aeroelastic analyses of a rectangular cantilever plate of varying aspect ratio is presented. The classical plate theory has been selected as the structural model. The main point that distinguishes this study from previously reported research is employing Peters’ theory to model aerodynamic effect which is not straightforward. The Peters’ aerodynamic model was originally developed to provide lift and moment, which is only applicable to the structural model based on the beam theories. In this study, using the basic concept of the Peters’ aerodynamic model in addition to utilizing the Fourier series, the pressure distribution is derived, which makes Peters’ model applicable to structural models based on plate theory. This combination provides a much simpler state–space aeroelastic model for plates in comparison to the prevalent panel methods, which could lead to a significant reduction in computational time. In addition, the aeroelastic response of the plate with respect to changes in the structural model from the beam theory to the plate theory is evaluated. By using data from an experiment carried out at Duke University, the theoretical results are evaluated. Furthermore, the differences in structural models obtained from the plate and beam theories can be divided into two distinct parts, which are responsible for differences in bending and torsional behaviors of the structure, separately. This approach enables us to measure the effects of differences of each behavior separately, which could provide with a new insight into the problem. It has been determined that the flutter speeds obtained from the beam and plate aeroelastic models are little affected by the difference in bending behavior, but rather is mainly caused by the difference in torsional frequencies.     

基于贝叶斯网络的船舶港口国监督检查滞留风险研究

针对港口国监督(Port State Control, PSC)检查的复杂性和不确定性, 基于贝叶斯网络理论构建船舶PSC检查滞留风险分析模型. 以东京备忘录(Tokyo MOU)中2014~2017年船舶PSC检查样本数据为基础, 运用R语言bnlearn包进行贝叶斯网络的结构及参数学习. 同时分别执行贝叶斯网络的正向、逆向推理, 定量表示各风险因素与滞留结果之间的相互作用关系, 找出导致船舶滞留的高风险因素, 实现不确定环境下船舶PSC检查滞留风险的全面动态分析. 实证表明, 模型具有较高的精确度, 可为检查人员的滞留决策及航运公司的安全风险管理提供有效依据.     

Measured and Predicted Solubility Phase Diagrams of Quaternary Systems LiBr-NaBr-MgBr2-H2O and LiBr-KBr-MgBr2-H2O at 298.15 K

The stable phase equilibria of quaternary systems LiBr-NaBr-MgBr₂-H₂O and LiBr-KBr-MgBr₂-H₂O at 298.15 K were studied by both experimental measurement(isothermal solution saturation method) and theoretical prediction(Pitzer model). The solubilities of the saturated solution have been determined experimentally and two stable phase diagrams and relevant water diagrams of the two quaternary systems were obtained. Results show that quaternary system LiBr-NaBr-MgBr₂-H₂O is hydrate II type as NaBr and NaBr·2H₂O coexistence. Its phase diagram consists of only one invariant point, four univariant curves, and five crystallization fields. The quaternary system LiBr-KBr-MgBr₂-H₂O is a complex type as the double salt KBr·MgBr₂·6H₂O formed. In addition to this double salt, the three single salts LiBr·2H₂O, KBr and MgBr₂·6H₂O also crystallize. In this paper, the solubilities of phase equilibria in above quaternary systems were also calculated by the Pitzer's electrolyte solution model. All the needed parameters can be obtained from the literature or be fitted by experimental data. On the Basis of the experimental and calculated results, the phase diagram of the quaternary system was plotted for comparison. It shows that the calculation results are consistent with the experimental ones.